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Abstract 

We consider the problem of integrability of the Poisson equations describing spatial mo- 
tion of a rigid body in the classical nonholonomic Suslov problem. We obtain necessary 
conditions for their solutions to be meromorphic and show that under some further restric- 
tions these conditions are also sufficient. The latter lead to a family of explicit meromorphic 
solutions, which correspond to rather special motions of the body in space. We also give 
explicit extra polynomial integrals in this case. 

In the more general case (but under one restriction), the Poisson equations are trans- 
formed into a generalized third order hypergeometric equation. A study of its monodromy 
group allows us also to calculate the "scattering" angle: the angle between the axes of limit 
permanent rotations of the body in space. 
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1 Introduction 



In some cases of the rigid body dynamics, in particular, in the problem of motion of a solid 
about a fixed point, the Euler equations for the angular velocity vector w € M'^ separate and 
can be integrated. Then, given a generic solution Lu{t), to determine the motion of the solid 
in space it is necessary to solve the reconstruction problem, that is, to find 3 independent 
solutions of the linear Poisson equations 

7 = 7X^(0, (1.1) 

7 g R'^ being a unit vector fixed in space. 

The most known example of solvable Poisson equations gives the Euler top problem, when 
generic ult) are elliptic (i.e., doubly periodic) functions and one particular solution "/{t) is 
also elliptic, whereas the other two are quasiperiodic (see, e.g., [71 [TO]). 

Similar, but formally more complicated solutions 7(i) appear in the case of the Zhukovsky- 
Volterra gyrostat (see p^[22] ). 

A nontrivial integrable generalization of the Euler-Poisson equations was found in [3], 
where the Euler equations have the standard form and (|l.ip are replaced by the equations 

7 — X uj(t), 

>t being an arbitrary odd integer number. It was shown that, like in the Euler top problem, 
the latter equations possess an extra algebraic integral, however a complete solution for 7 is 
still unknown. 

In present paper, following Suslov [TB;, we consider the motion of the rigid body about 
a fixed point in presence of constraint (w,a) = 0, a being a fixed vector in the body frame. 
Let I : M"^ R'^ be the symmetric inertia tensor of the body. The Euler equations for the 
angular velocity vector lo separate and take the following simple form 

^(Iw) = X + Aa, (1.2) 

where x denotes the vector product in R'^ and A is the Lagrange multiplier. Differentiating 
the constraint, we find A ~ — (Itj x uj, I^^a) / (a, I^^a). Therefore, p.2p can be represented as 

at (a, 11 ^a) 

which, in view of (w, a) — 0, is equivalent to 

^(Iw) = (Iw,a) w X I-^a. (1.3) 

In the sequel, without loss of generality, we assume that a — (0, 0, 1)"^, which, in view of the 
constraint, implies Lo-i = 0. This simplifies the Poisson equations to the form 

71 = -a;2(t)73, 72 =t^i(i)73, 73=^^2(^)71-^1(^)72- (1-4) 
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We also assume that the tensor I is disbalanced, i.e., is not diagonal in the chosen frame. (If 
a is an eigenvector of I, then all the solutions of (|1.3p are equilibria.) 

It is known that under these assumptions the system (jl.3p restricted to the plane UJ3 — 
has a line of equilibria points /13W1 + 123(^2 — and that the trajectories u>{t) are elliptic 
arcs that form the heteroclinic connection between the asymptotically unstable and stable 
equihbria (see Fig. 1 below). 

Whereas the reduced system (|1.3p is elementary integrable in terms of hyperbolic func- 
tions, it is believed that the corresponding Poisson equations (II. 4p are not, although we did 
not find a proof of that in the literature. A study of complex solutions of these equations 
was done in [TT]. A qualitative analysis of the behavior of 7(t) in the classical Suslov prob- 
lem, as well as in its multi-dimensional generalization was made in [21j , whereas some other 
interesting generalizations of the problem were studied in [8]. 

Contents of the paper. In Section 2 we present generic solutions of the Euler equations 
for the Suslov problem and formulate the problem of integrability of the system (|1.3p . (|1.4p 
by the Euler-Jacobi theorem. 

In Section 3 necessary and sufficient conditions of meromorphicity of solutions are ob- 
tained, they both require that one of the components /i3,/23 of the inertia tensor must be 
zero. 

Section 4 discusses general properties of solutions of the Poisson equations in the specific 
case /i3 = 0, which are compared with generic solutions of another famous nonholonomic 
system, the Chaplygin sleigh. It is also shown that the meromorphicity conditions on I are 
compartible with the restrictions on the inertia tensor of a physical rigid body. 

In Section 5 we show that in the case /13 = the Poisson equations can be transformed 
to a generalized third order hypergeometric equation. Using its monodromy group, we solve 
the classical problem of calculating the angle between the axes of limit permanent rotations 
of the body in space. 

When the sufficient conditions of meromophicity are satisfied, we observe that the corre- 
sponding hypergeometric series solutions reduce to products of polynomials and exponents, 
which are explicitly calculated. 

Next, we apply the differential Galois analysis to the Poisson equations and prove that 
when the parameters of the problem satisfy the necessary conditions of meromophisity, but 
not the sufficient ones, these equations and, therefore, the whole Suslov system, are not 
solvable in the class of Liouvillian functions. 

In Section 6 we present all the meromorphic solutions of the problem and the correspond- 
ing extra polynomial integrals in the explicit form. 

In Conclusion some relevant open problems are brieffy discussed. 

2 Generic solutions of the Euler and the Poisson equa- 
tions 

In the general case the components /i3,/i2 of the inertia tensor I are not zero, but, by an 
appropriate choice of the frame that preserves the constraint one can always make I12 = 0. 
For We also impose the normalisation det 1=1. Condition that I is positively defined means 
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that all main minors are greater than zero and we obtain In > and /11/22 > that gives 
I22 > 0. 



Then the Euler-Poisson equations (|1.3p , (|1.4p have the form 

oJi = hiihst^i + ^23^2)^2, = -hiihsuJi + 123^2)^1, (2.1) 

and 

71 = -W2(t)73, 

72=c^i(t)73, (2.2) 
73 = t^2(i)7i - ^^i(i)72. 
They have two first integrals, the energy and the trivial geometric one: 

Fi = Iniul + h2Lol. F2 = (7, 7> = 7i + 72 + 73'- (2-3) 



In this paper the main problem we consider is: For which values of lij the system (|2.1[) . 
(j2.2|) is integrable ? Here the integrability will be understood in the context of the classical 
Euler-Jacobi theorem, which relies upon the existence of an invariant measure and sufficient 
number of independent integrals. 

Although the system (|2.ip . (|2.2p does not have an invariant measure in the strict sense 
(the density of the volume form tends to as i — s- ±00), we can consider its restriction on 
the energy level Fi = E > 0, which consists of two open components. On each of them we 
choose time t and 7 as coordinates and then the restricted equations read 

71 =-^2(^)73, ^2 4) 

72 = i^l{t)j3, 

73 = ^^2(071 - ^i(i)72- 

They have the trivial integral F2, as well as the invariant volume form 

^ = dt A d7i A d72 A d73 , 

Thus, for the integrability in the Jacobi sense only one additional first integral F{t,'^) is 
required. 

The existence of such an integral is closely related to the existence of single- valued solu- 
tions of the Poisson equations. Namely, let P{t) — (Pi, P2, P3) a single- valued vector- function 
satisfying P ^ P x cu. Then the system possesses the single- valued integral 

F=(P(i),7), (2.5) 

Indeed, 

F = (P(<), 7) + {P{t), 7) = (P(0, 7> + {P{t),l XLu) = {P{t) +UJX P{t), 7) - . 

If, moreover, the components of P{t) are single- valued functions of the solutions a;i(<), 0^2(^)1 
then the system (|2.ip . (|2.2p admits an additional first integral 

F3(u;i,W2,7i,72,73) = (^'('^),7> 
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functionally independent with Fi and F2. 

We start with generic solution of the dynamic equations (|2.ip . which have the form 



At 



+ e- 



5 = 



13 
I 23 



(2.6) 



where 



AL 



23 



ci = ±2- 



13 



C2 = ±2- 



23 




^1^3-^22 + -^23^11 ' " ^1^3-^22 + -^23^11 V Al ' ' -^13-^22 + ^23-^11 

A being an arbitrary positive constant related to the energy integral. Note that for t — > ±00 
these expressions give points on the equilibria line /13W1 + /231-1J2 = 0, as required. 

Remark 2.1 For each A fixed, the Poisson equations (12. 2|) and the above functions Lo{t) are 
invariant with respect to time rescaling t t/A, which reduces the solution of to the 
form 



wiW = 

-^23 



e 



Cl 



^3-^22 + ^3^11 



e* + e * 

Cl 



C2 



±2/1 



e* + e 



C2 



6= - 
±2/23 



i23 



-^13^22 + ^23-^11 V -^22 



(2.7) 
(2.8) 



This imphes that, without loss of generality, one can study solutions of p.2p with the 
coefficients llJi,lu2 having the form (|2.7p . (|2.8p . This wi7i be assumed in the sequel. 

Note that choosing uj{t) in the form p.Tp implies that the energy integral is fixed to he 



Fi 



^13^22 



I 23 1 11 



(2.9) 



3 Painleve property 

In this section we investigate the following problem. 

Problem 3.1 For which values of the parameters lij all solutions of the system l|2.2p are 
meromorphic in C. 

The answer is contained in the following theorem. 
Theorem 3.2. All solutions of (|2.2p are meromorphic solutions iff either 

111 — I22 



/i3 ~ 0, and 
or I23 = 0, and 
where p is a nonzero integer. 



Illl22{Illl22l33 — 1) 

111 — I22 
-^1 1-^22 (-^11 ^22 -^33 — 1) 



(3.1) 
(3.2) 
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Figure 1: Generic trajectory u>{t) on the plane {u!i,u!2)- Black and white dots indicate stable and unstable 
equilibria. 



Below we also show that the second case p.2p can be obtain from p.ip by the linear 
transformation 

Wi ^ -W2, W2 ^ -Wi, 71^ 72 , 72 ^71, 73^73, ^gg^ 

III I22, I22 hi, hd. I23, I23 hs, I33 hi- 

Note that the conditions (|3.ip expressed in terms of the parameters (a, b, ci, C2) take the 
following form 

The proof of Theorem 13.21 is given in two next subsections. 



ci = 6 = and cj-Aa^^Ap^, p e Z* . (3.4) 



3.1 Necessary conditions — Kovalevskaya analysis 



Let us show that if all solutions of the system (|2.2p are single- valued, then one of the condi- 
tions (|3.ip . p.2p is satisfied. 

First observe that the right hand sides of the system (|2.ip . (|2.2p are homogeneous poly- 
nomials of degree 2 in the components of a; = (wi, 1^2, 7i, 72, 73)- Then we can apply the 
Kovalevskaya-Lyapunov analysis (see [TD] for details) to find that the system admits two 
solutions of the form 

x{t) = jd, where d = (cji,W2,7ii72,73)- (3.5) 
The first of them is given by 

1 _ 1 

^1 = T — ^ TT ;^ ^ ; ^2 = 



hihs + ^hsV I11I22 h3h2 — ^hsV hih2 (3-6) 

71 = 72 = 73 = 0, 
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and the second one is its complex conjugate. According to the Lyapunov theorem, if all 
solutions of (|2.2p are single- valued, then 

1. all eigenvalues of the Kovalevskaya matrix 

K{d)^^{d)+ld, 

are integers, and 

2. the Kovalevskaya matrix is semi-simple. 

For the solution d given by (|3.6p the characteristic polynomial of K{d) is 



(A-2)(A- 1)(A-H)1¥(A) 



det(if(d) - Aid) = ' 7 . .7 ; . . 2 , (3-7) 



where 

W{\) = (v/^/23 - i/l3 V^)[/22 + - /22(/l3 - i/^/23)'(A - 1) 

The roots of W{\) have the form 



, 1 I A A V (^11 ^ I22)W Illhs + 1-^13 V^22)'^ (VAl^23 ^ 1-^13 V-^22)^ 

Ai,2-1±A, A- ===——— — . 

V -'ll-'22(,-'llJ22-'33 ^ J^j 

If Ai and A2 are integer, then, in particular, they are real, which implies that 

/l3/23(/ll-/22)=0. (3.8) 

Moreover, if Ai, A2 G Z, then A'^ — , where p e Z*. This gives 

(ill - J22)Ull^23 ^ I22I13) ^ ^2 ^3 g-j 

.^11^22 (-^11^22^33 — 1)^ 

The condition p.8p gives three possibilities. Either /13 = 0, and this case gives (13. ip . or 
^23 = 0, which gives p.2p . or, finally, In = /22. In the last case, the Kovalevskaya matrix 
has the multiple eigenvalue A = 1 and is not semi-simple. As a result, we proved the 'if part 
of Theorem [ 



Remark 3.3 In terms parameters (a, b, ci, C2) conditions (|3.8p and (|3.9p iave the form 

aci+bc2 = and cj + 4 - 4:{a^ + b^) = Ap^ , peZ. (3.10) 



3.2 Sufficient conditions — analysis of the monodromy group 

The general solution (|2.7p of the Euler equations is single- valued. Thus, all solutions of the 
Euler-Poisson equations (|2.ip . (|2.2p arc single- valued if and only if all the solutions of the 
linear Poisson equations (|2.2p with a;i(<) and 0^2 (t) as in (|2.7p are single- valued. 
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Note that the only singular points of (|2.2p are simple poles located at to = ±f i niod iri. 
Then a necessary condition for the single-valuedness is that the monodromy matrices at all 
the singular points are identities. 

To show this, rewrite the system (|2.2p in the form 



7 = Mth, 



Ait) 



and observe that 
1 



t-to 



Aa + 0{t - to), Ao = 











-UJlit) 



-UJ2it) 
LOlit) 





(3.11) 



-b- 
a - 



|iC2' 

iici 



b — iic2 —a + iici 







(3.12) 



The monodromy matrix of the canonical loop around to is Mtg exp27riylo. A direct 
calculation shows that the eigenvalues of Aq are 



1 



pi = and P2.3 ^±^\Jcl+cl~ A{a? + b"^) + 4i(aci + feca). 



(3.13) 



Hence, if the conditions (|3.10p are satisfied, then A^ has eigenvalues pi = 0, p2 = and 
P3 — —p with p G Z*. In this case Mf^ = Id provided that there is no local logarithmic 
solutions in a neighborhood of tg. In order to check this, we use the following theorem (see, 
e.g., [U for the details). 

Theorem 3.4. Assume that in the linear system 



Y = B{t)Y, B{t) = t-^ ^ B,f 



(3.14) 



i=0 



the matrix coefficient Bq does not have a pair of eigenvalues such that their difference is a 
non-zero integer. Then there exists a matrix T given by a convergent power series 



1=0 



such that the linear map Y = TZ transforms (j3.14p into the following fo 



Z' = -BoZ. 
t 



(3.15) 



The form of the fundamental matrix Z Bo of (|3.15p is well known. Namely, let P be the 
similarity transformation reducing Bq into its Jordan form J — P~^BqP, where P is an 
invertible matrix with constant coefficients, and let Z j denote the new fundamental matrix. 
Then Z Ba = PZjP~^, and Zj has a block-diagonal structure such that the Jordan block 
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(of dimension v) corresponding to an eigenvalue A has the form 



1 ln< 
1 

.0 ... 



(^-1)! 







Int 
1 

Now assume that conditions p.ip or (|3.2p are satisfied. By the existence of the hnear 
transformation (|3.3p , without loss of the generality we can assume that conditions (13. ip are 
satisfied. 

Expressed in terms of parameters (a, 6, ci,C2), the latter take the form (|3.4p . or, setting 
C2 = — 2co, = /11//22, the following form 

(c2-l)a2=p2. (3.16) 

We also denote (fi — — I. 

Now we check the presence of logarithmic terms in the solutions. 

Theorem 3.5. The Poisson equations (|2.2p with uji given by (12. 7p with Ci = 6 = 0, C2 = 
—2ca = — 2V (P + Ip/d, p ^ li, do not have logarithmic terms in local solutions around 
singular points to ■ 

Proof. For ci = 6 = 0, C2 = —2ca = + Ip/d, the matrix Aq reads 



ipVd^ + 1 
d 



\p\/ + 1 p 



(3.17) 



Using the similarity transformation S , one obtains its diagonal form: ^ S ^ A^S with 






















0" 




d 


io = 











, s = 


1 














d 



1 



Vd^TT 
1 




iVd^TT 
1 

1 



Thus we see that all the eigenvalues of Aq are integer and their differences are also nonzero 
integers. Under this transformation the matrix A of the linear Poisson equations takes the 
form 

ipsech {t){d'^ + 1 + i sinh(i)) ip(sech (t) + i tanh(i)) 



A 



{d^ + l)p(-i sech {t) + tanh(i)) 






(d^ + l)p(-i sech {t) + tanh(i)) 
d3 ■ 
ip(sech (t) + i tanh(i)) p[-i{d^ + 1) sech {t) + tanh(i)] 



2d 




2d 



d? 
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Transformation 7 



r defined by 7 = UT with 

2p 



U ■ 



ITT 
t 

2 



ITT 
t 

2 





gives tlie equivalent linear system T = B{t)T with 

B{t) = U^^AU -U^^ 
This matrix has residual part diagonal 



U. 





f-p 








Bo = 












^ 





-p 



which satisfies the assumptions of Theorem 13.41 Since all the Jordan blocks of Bq are one- 
dimensional, the logarithmic terms do not appear in the solutions of the Poisson equations. 

□ 



4 Special case of the inertia tensor. General properties 
of solutions 

Assume now that, apart from I12 = 0, the first meromorphisity condition in (|3.ip is also 
satisfied, that is, I13 = 0, which imposes an essential restriction on I. 
Then the reparametrized solutions (|2.7p . (|2.8p are reduced to 

a (e* — e~*) , , 2ac , , , 

e'^ + e ^ e'^ + e 

1211?,?, — I^? fl- 



22-'33 — ^23 /^ll 
-'23 V -'22 

It is seen that they have equilibria along the line uj2 — 0. 

Using the parameter p introduced in (|3.ip . and setting d — p/a, we can rewrite the whole 
system ((2T]) . ([2?2)) in the form 

d 2 
p(d'' + 1) 
d 

UJ2 = ^1^2, 

P (4.2) 

71 = -Ci^273, 

72 = W173, 

73 = <^2l\ - ^ll2, 
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which always has two first integrals 

^ [d^ + 1)loI + ul , F2 = 7? + 72 +1I (4.3) 

Note that on the solutions (|4.ip one has Fi = a^c? = As follows from Theorem 13.21 the 
solution of this system is meromorphic if and only if 




^11^22 (Al-^22-?^33 — 1) d 



is a non-zero integer. 

Since d — \/(? — \ and p can take positive, as well as negative values, we choose sign + 
in the first equation in (|4.2p . 

Proposition 4.1. For any p > there exists a "physical" rigid body with the inertia tensor 
I and /i3 — such that its eigenvalues satisfy the triangular inequalities and its components 
satisfy the meromophisity condition H3.1\) . 

Proof. In view of the structure of the tensor I, for its eigenvalues to be positive, the diagonal 
entries In, I22, ^33, as well as I22I33 — lis^ must be positive. 

Assume that In > I22 > and 133 + I22 > In- Then, given such /ii,/22,^33, the 
condition p.ip represents a quadratic equation with respect to /I3, whose solutions are 

4 = 2{Iu~l22) ^ " ^''^ ^ pVp'+4/33(/i1-/22)) • 

One can check that both solutions are positive, hence the obtained I23 are real. Let us chose 
the solution with sign — at the square root. Then we get 

/22/33 - /I3 = ... , ^ (V+Pv/p'+2/33(/ll-/22)) . 

^i^ll — J22j ^ ' 

Since p and In — I22 are positive, we have I22I33 — I23 > 0, as required. 
Next, the eigenvalues of I are Iii,J2,J3, where J2, J3 are the solutions of 

- (^22 + /33)y + ^22^33 " I23 = 0- 

The determinant of the latter equation is D — {I22 — 133)^ + 4/I3 > 0, hence both roots are 
positive. 

Finding the roots, one can also check that if^ — ( J2 — J3)^ > 0, which gives the triangular 
inequalities In + J2 > J3, hi + J3 > J2- The third inequality follows from the assumption 
I33 + I22 > hi and the fact that J2 + J3 = h3 + l22- □ 

Steady-state rotations of the body in space. Comparison with the Chap- 
lygin sleigh. Solutions of the Poisson equations (|1.4p with the coefficients u!{t) in (|2.7p 
describe the evolution of the rigid body in space, which is an asymptotic evolution from one 
steady-state rotation to another one. The equilibria of uj{t) correspond to the steady-state 
rotations themselves. 
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Figure 2: A generic trajectory 7(t) — {gl, g2, g3)^ with the boundary condition (|4.4p . 



In the particular case /13 = 0, as i — > ±00, the motion in space tends to rotations about 
the axis (1,0,0) fixed in the body with the angular velocities =Fa. Note that the spatial 
orientations of the axes of the above steady-state rotations are generally not the same. 

This can be illustrated by solutions of (|1.4[) with the initial (or, rather, boundary) con- 
dition 

lim 7(t) = (-1,0,0)^. (4.4) 
t — >— 00 

As follows from the structure of (|1.4p and (|4.ip . in this case, as i — > -l-oo, 7(t) tends to the 
periodic trajectory 

(cos(Ai/'), sin(A?/') cos(at), sin(Ai/') sin(at))"^, 

where Atjj is the angle between the above axes of the steady-state rotations in space. That 
is, the trajectory 'y{t) in the body frame tends to a uniform rotation about the axis (1, 0, 0) 
(see an example in Fig. 2). 

A similar behavior occurs in the "non-compact" version of the Suslov problem: the 
Chaplygin sleigh, a rigid body moving on a horizontal plane supported at three points, 
two of which slide freely without friction, while the third is a knife edge (a blade), which 
allows no motion orthogonal to its direction (see for the details, among other, [11 [12]). The 
configuration space of this dynamical system is SE{2), the group of EucHdean motions of 
the two-dimensional plane R^, which can be parameterized by the angular orientation 9 of 
the blade, and the position (x, y) of the contact point of the blade. The presence of the blade 
determines a nonholonomic distribution on the phase space T SE{2). 

Like in the Suslov problem, the equations for the two linear momenta and the angular 
momentum in the body frame separate and are solvable in terms of hyperbolic functions. 
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In contrast to what is known about the solvabihty of the Poisson equations (|1.4p . the 
kinematic equations for the Chaplygin sleigh are integrable for any initial conditions [12]. 
In particular, as t —^ ±(X), the motion of the sleigh tends to a straight line uniform motion, 
however along different directions. (This can be compared with the steady-state rotations 
of the body about different axes in the Suslov problem.) A typical trajectory of the contact 
point of the blade is shown in Fig. 3. In addition, in (12j the following remarkable property 
was observed: the angle between the limit straight line motions of the sleigh does not depend 
on the initial conditions, but only on dynamical parameters of the sleigh. 

In this connection the following classical questions arise: is it true that the angle Aip 
between the axes of the limit steady-state rotations of the Suslov problem depends only on the 
components of the inertia tensor I (that is, it does not depend on the energy of the motion) 
? And, if yes, how Aip depends on the moments of inertia lij ? 

It appears that the answer is positive and the function A'ip{Iij) can be calculated explicitly 
without solving the Poisson equations (see Theorem 15.41 in Section 5). 




Figure 3: A generic trajectory of the contact point of the blade on the plane {x,y). 

5 Solvability analysis 

5.1 Solvability in the class of generalized hypergeometric functions 

In this subsection we again assume that /13 — and fix the energy the level to (|2.9p . Hence, 
LUi{t) and uJ2{t) are given by (|4.ip . with 

a = ^, d^=c^-l. (5.1) 
d 

Note that here we do not assume that p G Z*. 
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The Poisson equations (j2.2p can be rewritten as the foUowing third order equation 

ao'7i + ai7i + 0271 + a37i = 0, (5.2) 

where 

ao(i) = {(f + ai{t) = dp{2{(f + 1)ujI - lu^), 



a2{t) = {<f + l)uji[{d? +P^)loI +P^ujI], asit) = -dpujl[{d? + 1)ujI + ujI 
Now introduce new independent variable 

, 4 



(5.3) 



Then 



. _ 2dzuji .. _ 2z(2p2 + 2dV - 3(1 + d^)p^z) 
^~ ~' (d2 + l)p2 ' 

... _ 8dz(p2 + _ 3(-^2 _^ i)p2^)t^^ 

and (15. 2p . ()5.3p transform to 

7^ + &i7i + b2l[ + ball = , (5.5) 
where the prime denotes the differentiations with respect to z, and 



u 2 1 ^ (d2 + l)p2 



z z-1' 8d2(z-l)2z2' 

p2(z - 1) + ^2(1 + z{-6 - p^{z - 1) + 4z)) 

4^2(7^1)2^2 



(5.6) 



Now, setting 71 — \/ z ~ lu(z), we obtain the foUowing equation for the function u{z) 



z2(l-0)u"' + i(4-9z)zw" + 



1 



w' + -(p2 - l)u = 0. (5.7) 



This is a generahzed third order hypergeometric equation, whose canonical form is (see e.g., 

m) 

z^{\ - z)u"' + z[l + /3i + /32 - (3 + ai + a2 + a3)z]u" 
+ — (1 + ai + ^2 + 03 + aiQ;2 + a2Q;3 + a3ai)z]u' — aia20tau = . 

One can identify (|5.8p and (|5.7p by setting 

1 1+p 1-p d-\p d + ip , . 
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Equation (|5.8p belongs to the class of generalized hypergeometric equations of order n for 
n = 3. If we denote by 6* = zd/dz, then this class takes the form 

[(0 + /3i - 1) • • • (0 + Pn-i - 1)0 - z{6 + ai)---{e + a„)]u = 0, (5.10) 

see e.g. [Ul [2]. Recall that one of its solution is the generalized hypergeometric function 
introduced by Thomae [17] 

p (ai, an \ {ai)k ■ ■ ■ {an)kz'' 

Here (a)^ — a{a + 1) • • • (a + fc) is the Pochhammer symbol and (a)o — 1. 

It is known (see e.g., |13j) that equation (|5.8[) has three regular singularities over CP^ 
with the exponents 

0,l-/3i, 1-/32, at z = 0, 

0, 1, /3i + /32 — ai — Q!2 — as, at z = 1, 
Q!i,a2,Q;3, at z = oo. 



and if /3i , /?2 , /?i — /32 ^ then the general solution of (|5.7|) around the origin is given by 
linear combination 

u = CiTi{z) + C2T2{z) + C3.F3(z), (5.12) 

^l(^) - 3i^2 (^^^ , 

■77/ \ i-/3i i;, /ai-/3i + l a2-/3i + l as - /3i + 1 i /r;iQ\ 
^2(z) = z ^^^3^2^ 2_^^ ;z) , (5.13) 



_ p / ai - /32 + 1 a2 - /32 + 1 as - /32 + 1 

' /3i - /32 + 1 2-/32 



Ci, C2, and C3 being constants of integration. 

An immediate observation is that if one of in (|5.1ip is a negative integer, then the 
series truncates. In view of (j5.4p . this happens precisely when p is an odd integer, and then 
the particular solution J-i{z) takes the form 

TAz) = 1 + ' (5 14) 

j^, (2j)!! (d^+p^)---m-i)d^+p^)' ■ ^"-'^^ 



Since 71 = ^/]~~zu{z) and z = 4/(e* + e^*)^, this corresponds to a solution 71 (i), rational 
in e* and meromorphic in having poles of order p at f = 7ri/2 (mod 7ri), as predicted by 
Theorem [ 



Transformations for ^2 (^) 5-^3 (-2) • Due to the symmetry in the definition (|5.9p of 
the parameters a,/3, the above series !F2{z) ^ !F'i{z) can be expressed in terms of customary 



15 



hypergeometric functions 

F (a, ^; = 1 + + 
of the variable y — e^*. 

Proposition 5.1. Under the condition \5. 9i) one has 

dzLip 

(i + y)p-i(y-i)^' 

where 

2d 



•^2,3(2) = -— — r—r-, -ry 2d ■ F{a±,b±; c±; ~y)F{a±,b±; c±; ~y), (5.15) 



z = 4y/(l + 2/)^ >f± 



(d±ip)(3d±ip)' 



F = 1)?/- dTip]-F(a±,6±;c±;-2/) - 2dy{y + l)-^F{a±,b±;c±,-y) 

and 

Proof. In view of (|5.9p . the last two formal solutions in (|5.13p can be written as 
.^2,3(2(2/)) = (4y) (y+l) 3F2(^„ .^^"^^ , 1 ' ■ 



2a± + 2/3± - 1 a± + /3± + i ' (y + 1)2 y 

The functions 3F2(') of such kind split into a product of two customary hypergeometric 
functions 2^i(') = F{'), namely 

2a 2/3 a + 13 

1-2/3-1 a + /3+i 

where 7 = a + /3 + 5, see e.g. formula (8) on page 86 in _6J. In our case 



I 2a + 2/3-1 a + /3+i ] = F {a, P;r,x) F {a, (3;-/ - l;x) , (5.18) 



a = a±, ^ = /3±, 7± = a± +/3± + i = ^ ± a; = 4y/(7/ + 1)^. 

Next, due to the special form of 7, the first hypergeometric function _F (a, /3; 7; x) can be 
written as series in y by using the following quadratic transformation 

a a 1 , ,4?/ 

-,-H bA+a-b,- 

2'2 2 ' (1 + y)^ 



F[-,- + --b,l + a-b,j^^] ^il + yrF{a,b,l + a~b,~y), (5.19) 



see e.g. page 64 in [6]. Now identifying a± — a±/2 and /3± = a±/2 + 1/2 — fe±, we obtain 

2-p,ip , 1-p , 3 ip 

a+ — ± — , 0+ — , c+ — I + a+ ~ 0+ — — ± — , 

^ 2 2d' 2 ' * =^ =^ 2 2d' 
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which coincide with (|5.16p . Inserting these parameters into (I5.19p . we get 

F(a±,(3±;j±;^^j^) = {y + 1^^ F{a±,b±; -y). (5.20) 



It remains to apply a similar quadratic transformation to the factor F{a±, /3±; 7± — 1; 4^//(J/+ 
l)■^). To do this we first use the relation 

d" 

(7 - n)„ z-'-'-^Fia, (3-j-n;z)^ — [z-'-'Fia, p- 7; z)] , 

(7 — n)n being the corresponding Pochhammer symbol (see e.g. equation (22) on page 102 
in [B]). Setting here n = 1, we obtain 



1 z d 

Fia±,f3±;-f± - l;z) = — F{a±, l3±;-f±; z) H -T:^F{a±, (3±;j±; z), 

7± 7±(7±~l)dz 

4y 



(5.21) 



(j/ + l)2- 
In view of the chain rule 



(y+l)3 d 



dz 4(y-l)dy' 
and (|5.17p . the expression (|5.2ip takes the form 

4y \ 2d{y + 1)*^* 



F a±,/3±;7± - 1; ■ 



{y + l)y {d±ip)i3d±ip)iy~l) 



[d{p- l)y - dTip]F{a±,b±;c±,-y) ~ 2dy{y + l)-^F{a±,b±;c±,-y) 

dy 

Substituting this, as well as (|5.20p . into (|5.18p and simplifying, we finally get the expres- 
sion (HjH). □ 

Note that whenp is odd positive integer, the parameters b± — (1— p)/2 become negative 
integer, and the hypergeometric series F{a±,b±; c±, —y) in ()5.15p convert into polynomials 
of degree {p — l)/2, 

F(a±, c±, -y) = ''f:' ^"fft^' {-vf ■ (5-22) 

Next, after a simplification, the series J^(a±, fo±; c±, — y) in (|5.15p becomes the following 
polynomial of degree [p — l)/2 (the term of the highest degree (p — l)/2 + 1 annihilates): 

p. , N ^'^y-J'^ {a±)k{b±)k {2dk + d + pi){Adk~ dp±ip) 
na±,6„c,,-y)= g dp-2dkTip 
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As result, we conclude that in this case, 

^My)) = (i + ^)"1(^_,) y'^P^iy), (5-24) 

where P±{y) — F{a±,b±;c±,~y) F{a±,b±;c±,—y) are polynomials of degree p — 1 with 
complex conjugated coefficients. 

The fact that in this case the series JFi is a rational function of y follows from the formula 
(I5l4l) . 

The explicit form of the corresponding mcromorphic solutions of the Poisson equations 
(|2.2|) for odd p is given in Section \6\ 



Monodromy for J^i (z) and the angle Aip. Now we show how using the monodromy 
group of the equation (|5.10p , which was analized in [Kl [11 [S] , one can calculate the angle 
Alp between the axes of the limit steady-state rotations of the body in space for generic 
p e M. As follows from (15. 4|) . the values t = ±oo, t = 0, and t — 7ri/2 (mod ni) correspond 
to the singular points z = 0, z = 1, and z = oo respectively. As time t evolves along the 
real axis from — oo to oo passing by t = 0, the variable z makes a loop in C embracing 1 in 
positive or negative direction. 
Next, we need 

Lemma 5.2. Under the substitution (j5.4p . the component 71 of the special solution j{t) G S*^ 
satisfying the boundary condition (j4.4p has the following form around the origin z — 

^i{t)^-^l-zit)Ti{zit)) 

- ./r^ F ( (l+P)/2 {l-p)/2 \ 

= -VI - ^ 3F2 (^^^2 - ^p/{2d) 1/2 + ^p|{2d) ' ' ^^'^^^ 

Proof. As follows from (|5.12p . the most general possible expression for the component 71, 
which ensures limf^_oo7i — —1, is 

71 (0 = -VI -Z {Tl{z) + C2T2{Z) + C3.F3(Z)) , 

C2,C3 being arbitrary constants. On the other hand, from the Poisson equations and the 
solutions (|4.1[) one has 

71 4 ^* ~ '^'^^ — ^* ~ '^'^^ 

LOi ac{e*' + e^*)^ rfz ac dz 

The above special solution requires that limt^_oo73 = 0. Then, setting here 71 =71, 
differentiating J'i^2,d,{z) and taking into account the values of /3i, /32, we see that this condition 
holds if and only if C2, C3 above are zero. □ 

When z makes a loop around 1, the root VI — z changes sign, whereas the hypergeo- 
metric series J^i{z) transforms according to the following proposition, which is a corollary of 
Theorem 4.1 in pj]. 
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Proposition 5.3. When z e C makes a loop around the singular point z — 1 in positive 
direction, the solution J'i{z) around the origin undergoes the monodromy 

Ti- 2i exp(^i + /32 - ai - 0^2 - as) {<y\T\ + 0-2-^2 + 0-3-^3) (5.26) 
sin(7rQ!i) sin(7rQ;2) sin(7ra3) 



0\ 



0-3 = 



sin(7r/3i) sin(7r/32) ' 
sin(7r(/3i - ai)) sin(7r(/3i - 02)) sin(7r(/3i - 013) 

sin(7r/3i) sin(7r(/3i - /32)) 
sin(7r(/32 - a.\)) sin(7r(/32 - a2)) sin(7r(/32 - as) 
sin(7r(/32 - sin(7r/32) 



Now combining formula (|5.25p with the proposition, taking into account (|5.9p . then eval- 
uating JFi, JF2, •?^3 at z = 0, one concludes that 

. „. / . ,„x . „ cos^(7r»/2) 
hm 71 = 1 - 2iexp 7ri/2 ai = 1 - 2 ^,,.Ji,'.y 
t^oo cos^((7rip)/(2q)) 

Since limt_^_oo 7i = ^1 and cos(AV') — — limt^oo 71, we obtain the following 

Theorem 5.4. The angle A?/; between the axes o] the limit steady-state rotations of the body 
does not depend on the energy and, when /13 — 0, is uniquely defined by relation 

A-ib cos(7rp/2) 
cos ■ — 



2 cos((7rip)/(2g)) 



where, as above. 



hi — 122 / hi 



I22 ' y A 1^22 (-^11-^22-^33 — 1) 

and it is assumed that dot I = 1. 

Note that when the parameter p is odd integer, Aip is always tt, regardless to value of q. 

We also add that the formula of Theorem 15.41 stands in a perfect correspondence with 
numerical integration tests. 

It remains to study the case of even integer p, when, according to Theorem 13.21 all the 
solutions of the Poisson equations (j2.2p are meromorphic, but are not given by truncated 
generalized hypergeometric series. This will be one of the subjects of next subsection. 

5.2 Differential Galois analysis 

Here we recall one result of G. Darboux which was formulated in Chapter II of his Theorie 
generale des surfaces. 

Lemma 5.5. Assume that ^^^"^ = ^^"^"^{t) is a real solution of the Poisson equation 

j = jxuj{t), (5.27) 

where Lj(t) is a real vector, and ^^^^ satisfies {'^'^^\^^^^) = 1. Then the remaining two 
solutions of ()5.27p linearly independent with 7'^"'^^ can be found with by a single quadrature. 
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Proof. Let us restrict equation (|5.27p to the unit sphere (7,7) = 1. As coordinates on it we 
choose 

73 + 1 7i + i72 
Ui = : — , U2 = — , (5.28) 

71 - 172 73 + 1 

so 

I-U1U2 .I+U1U2 U1+U2 /r- or,^ 

7i = , 72 = 1 , 73 = • (5.29) 

Ml - 1*2 Ul - U2 Ui - U2 

Now, it is easy to check that ui and U2 satisfy the following Riccati equation 

ii^A + Bu + Cu^, (5.30) 

where 

A= i(w2(t)-ic^i(i)), B = -ic^3(t), C ^]^{w2{t)+iwi{t)) (5.31) 

Formulae (|5.29p show that knowing two different solutions of (I5.30p we determine one solution 
of (|5.27p . On the other hand, having one real solution 7^^^ of (|5.27p we have two different 
solutions of Ricatti equation (|5.30p which are given by formulae (|5.28p . 

Let uq be a solution of (I5.30p . Then, as it is easy to check ui — — 1/uq, where z* denotes 
the complex conjugate of z, is also a solution of this equation. It is well known that the 
general solution u of a Ricatti is determined by its three different particular solutions uq, ui 
and U2- Namely, we have 

^u-^(u^-^^^^ (5.32) 

[U - Ml) [U2 - Uo) 

where Co is an arbitrary constant. From this formula we deduce that knowing only two 
different solutions uq and ui, the general solution can be obtain by a single quadrature 



u — Ul 



Co exp 



C(s)(uo(s) - Mi(s))ds 



(5.33) 
□ 



Putting u = —v/C transform equation (|5.30p to the form 



C 

V = Ai+ Biv-v"^, where Ai = -AC, Bi = B + —. (5.34) 
Then v = w/w where 

w~Biw-Aiw^O. (5.35) 

The general question is whether we can find an explicit form of solutions of Poisson 
equation (|5.27p for given functions uji(t). A proper setting to this question is given by the 
differential algebra. Namely, we assume that uji {t) are elements of certain differential field L 
with C as a subfield of constants. In the considered Suslov problem L = C(e*). 

Let if D L be the Picard-Vessiot extension for equation (|5.27p . We say that the equation 
is solvable iff extension K D L is a. Liouvillian extension. In this case all solution of the 
equation are Liouvillian. By the known Kolchin theorem, the extension K D Lis a. Liouvillian 
extension iff the identity component of the differential Galois group G{K/L) is solvable. 
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Proposition 5.6. Assume that equation (|5.27p has a solution Liouvillian over L = C(e*). 
Then all its solutions are Liouvillian. 

Proof. Let 7^^-' be a Liouvillian solution of (|5.27p . Assume that (7^^^ 7*'"'^'') = I. From the 
proof of Lemma 15.51 we know that to find other solutions of (|5.27p it is enough to find a 
general solution of the Ricatti equation (|5.30p . But 7*^^-' gives us one Liouvillian solution uq 
of equation (|5.30p . The general solution of (|5.30p is given by 

M = Mo H — , (5.36) 

y 

where y satisfies linear equation 

y= -{B + 2Cuo)y-C. (5.37) 

Hence y is Liouvillian and all solution of (|5.30p are Liouvillian, and thus all solutions of (|5.27p 
are Liouvillian. 

If 7'^-'^) is a Liouvillian solution of (I5.27P and (7^-'^^ , 7'"'^'') — 0, then, as uji{t) are real, Re7^"'^-' 
and Im7(^) are linearly independent real solutions of (|5.27p . Moreover 

7 = aiRe7(i) +02^x17(1) + a^iRe'^^^'^) x {Im-f^^'^) (5.38) 

is a solution of (|5.27p for arbitrary ai, a2 and as. In effect all the solutions off (|5.27p are 
Liouvillian. □ 

From the above proved fact we obtain the following. 

Proposition 5.7. Poisson equation (|5.27[) has a Liouvillian solution iff all solutions of 
linear equation (j5.35|l are Liouvillian. 

Proof. If equation (|5.27p has a Liouvillian solution, then, by Proposition l5.61 all its solutions 
are Liouvillian and formulae (|5.28p show that all solutions of the Riccati equation (|5.30p . 
as well as of the transformed Riccati equation (|5.34p . are Liouvillian. As v — w/w, w is 
Liouvillian, and thus all solutions of (|5.35p are Liouvillian. 

On the other hand, if equation (j5.35p has a non-zero Liouvillian solution, then all solutions 
of Riccati equation (|5.34p . as well as of (|5.30p are Liouvillian. Thus, by formulae (|5.29p . 
Poisson equation (|5.27p has a Liouvillian solution, so, by Proposition 15. 6[ all its solutions 
are Liouvillian. □ 

For the considered problem 7 = 7 x a;(t), equation (|5.35p reads 

w" +p(z)w' + q(z)w ^0, z = e\ '=^, (5.39) 

dz 

where 

+ 4icz - 4) - 1 



p{z) 



Z(z2 + l)(22 + 2icz-l)' 



p2 2^ p2 



(5.40) 



4(c2-l)z2 (^2 + 1) 
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This equation is derived under assumption that relation (j3.16p holds true. In what follows 
we will work with the reduced form of the above equation, i.e., with 

PI 1 

y" = r{z)y where r{z) = q = 2^'^^) + 4^^^)^ ~ (^'^-^^ 
where P is a polynomial of eight degree 

8 

P^Y.p,z\ (5.42) 

with the following coefficients 

PQ^Ps^d^ +p^, pl^pj^4ic{2d^+p^) , p2=pe^~4{M^+p^) , (5.43) 

pI - p5 - 4ic ((4p2 - 2) + p-^) , (5.44) 

P4 = -2((-8rf2 + 8(d2 + 2)p2 + i)rf2_^5p2j^ d^^c^-l, (5.45) 

and 

Q = 4 - 1) z2 (^2 ^ ^^2 ^2^^ _ i (^2 _ ^))2 ^^^^^^ 

Note, that we assumed that c 7^ ±1, i.e., taking into account relation p.l6p . p ^ 0. This 
case we consider later. 

From (|5.46p it follows that equation (|5.4ip has six regular singularities Si, i = 0, . . . 5, and 
sq — 0, Si — S2 — i, S3 = — i(c+ d), S4 = — i(c— d), and S5 = 00. The respective differences 
of exponents at these points are following 

Ao = A5-i^, Ai=A2=p, A3-A4 = 2. (5.47) 
d 

We observe immediately that at four singularities si, S2, S3 and S4 differences of exponents 
are integer thus at local expansions of solutions around of this points logarithmic terms 
can appear. Here we only sketch how this happens, for more details see e.g. [20]. For 
linear equations of the second order local solutions yi and 2/2 around a singular point Si are 
postulate in the form infinite series with leading term that are equal to the exponents at this 
point. 

In the case when differences of exponents A^ = pj^'' — /op"* is integer such two expansions are 
in general not functionally independent and then 7/2 is constructed in different way from yi 
by a quadrature that usually leads to logarithms. 

By direct calculations one can check that solutions around S3 and S4 do not have such 
terms independently on value of p. For singular points si and S2 direct calculations for small 
p show that the following conjecture is true 

Conjecture 1. Singular points si and S2 are non-logarithmic for any odd integer and loga- 
rithmic for any even p G Z* . 
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For singular points si and S2 the difference of exponents depends on p and when p grows 
we have calculate the expansions of wi with more and more terms. By this reason we 
cannot to check the presence of logarithmic terms for an arbitrary p but we can make such 
calculations effectively for any chosen value of p and we did this up to p = 10. 

The presence of logarithmic terms restricts very strongly the possible forms of elements 
from the differential Galois group, thus at first up to end of this section we will consider the 
case p even and we prove. 

Proposition 5.8. The differential Galois group of equation (|5.4ip for even p ^ Z,* is 
SL(2,C). 

Proof. The presence of logarithms means that the differential Galois group of the considered 
equation cannot be neither a subgroup of the infinite dihedral group (because it contains a 
non-diagonalizable element), nor a finite group. If it is contained in triangular group, then 
by the same reason it cannot be its proper subgroup i.e. diagonal subgroup. Thus wc have 
two possibilities that it is contained 

1. in the whole triangular group, 

2. is SL(2,C). 

If the first possibility occurs, then equation (|5.4H) has an exponential solution of the form 

4 

y ^ R{z)W{z ^ Si)P\ where i?(z) e C[z], (5.48) 

1=0 

and Pi £ {Pi^\ Pp"*} are exponents at points Si, for i = 0, . . . ,4. Expanding this solution at 
infinity we find that 

5 

n = — ^ Pi > 0, where n — degR. (5.49) 

i=0 

We have 

= i(l±A,) for z = 0,...,4, (5.50) 

and 

P5--^(1±A5). (5.51) 

The above implies that we have to choose po and p^ such that po + ~ 0. As point 
si and S2 are logarithmic, then the only choice is p^^^ = pj^^ = (1 +p)/2 for p > 0, or 

(2) (2) (1) (2) 

Pi P2 (1 ~ -P)/2 for P < 0- At S3 and S4 both exponents pl — 3/2 and p- = — 1/2 
for i = 3,4, are possible. Thus for p > 0, n G {—p,—p~ 2,— p — 4} and for p < we 
have n € {p,p — 2,p — A} but all these admissible values are negative. This shows that the 
equation does not have an exponential solution (|5.48p . Thus, the differential Galois group of 
the equation is SL(2,C). □ 

By the Kolchin theorem this means that in this case the Poisson equations are not solvable 
or, more precisely, the following theorem holds. 

Theorem 5.9. Euler-Poisson equations in the meromorphic case defined by (|3.4p for even 
p € 1,* are not solvable in the class of Liouvillian functions. 
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6 Explicit meromorphic solutions and first integrals 
6.1 Meromorphic solutions 

As was observed in subsection 15.11 for odd p € N the third order hypergeometric equa- 
tion (|5.7p for the variable u = ji/^/z — 1 has three independent quasi-polynomial solutions 
•^1,2,3(2) given by ((ETi)) and ((Oi)) . that is, 

4 



(e* + e-*)2' 



2/ = e , 



where Q{z), P±{y) are polynomials of degree (p — l)/2 and p — 1 respectively. 

Now, taking into account the relation between the solutions w{z) of this hypergeometric 
equation and those of the Poisson equation (j5.2p . we arrive at 

Theorem 6.1. For odd p e N the third order equation (|5.2p for 71 (t) has independent 
meromorphic solutions 

W = j^^, '"ji^kie' + e-)— , (6.1) 
^ (2fc-l)!! d^^(p^-l)---(p^-(2fc-l)^) 

(2fc)!! (d2+p2)...((2fc-l)d2+p2) ' "0 
^ ' k=0 ^ ' 

where (•)* denotes the complex conjugation, and the coefficients hi are uniquely determined 
from the polynomial product 

V^, k ^^V^^^ {a)k{h)k I ^fc ^^^^^ (a)fc(b)fc (2d/i; A- d + pi){Adk - dp + ip) 

2-p ip 1-p 3 ip 

a = , = , c ~ — . 

2 2d 2 ' 2 2d 

The proof is straightforward: Substituting the expression (|5.14p into 

yields (|6.ip . Next, in view of (|5.24p . the product 

e* - e-* 

7 = -T——[^^Aziy)) 

e^ A- e ^ 
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gives dU), (EJl). □ 

Note that all these solutions satisfy 

7['\~t) = l['\t), hm j['\t) = Tl, ^ lim 7f ''^ = 0. (6.4) 

t — '±00 t — >±oo 

Now let 7*^*-' = (Ti^'j 72*') 73*')"'" be the corresponding vector solutions of the Poisson 
equations (jl.4[) . Given as in Theorem 16. 1[ the two remaining components of can be 
calculated by differentiations, using the system (|1.4p . 

Namely, let us write the above meromorphic solutions in form 

'1 /-I I ™2^T)' '1 



(1+X2)P' '1 (1+X2)P' 

Pi(x), ^2(2;) being polynomials of degree 2p and 2(p — 1) respectively. Then, in view of p.4p . 
one has 

73 = , and then 72 = () = — . (6.5) 

UJ2 ^1 CLX 

Applying these formulas and using the expression (|4.ip for Ci;i,Ci;2, we get 
(1) ^ d{{\ + x-')P[{x) - 2pxP^{x)) 

2pVrf2TT(l+.T2)P 

(1) _ x[d2(x2 + l)((a:2 + l)P^{x) + (2 - ^p)xP[{x)) + 2p{2p + (P{2p - l){x'' + l))Pi{x)] 
'^^ ~ 2p2Vd2TT(l + 2;2)P(a;2-l) 

and 

(2) _ x-'-^[dx{x^ + 1)P^ + [d-ip+ [d-ip-2dp)x'^)P2] 
73 — 



2pVd2 + 1(1 + a;2)P 



(2) X dW 

72 



2p^Vd^Tl{l + x2)p(a;2 _ 1) ' 
W = -dx{x^ + l){dx{x^ + l)P2ix) + 2{d ~ 2d{p - l)a;2 _ \p(x^ + l))^'2(a;)) 
+ (p2(a;2 - 1)2 _ 2^2(1 +p(2p - 3))x2(x2 + 1) - idp{x^ + 1)(-1 + {4p - 3)x^))P2{x) 

(3) (2) 

Expressions for 7^ are complex conjugations of those for 7^ . 

It remains to mention that all the components of 7^*-' have poles of order p att = |-i+7riiV, 

N eZ. 



Remark. Note that for real moments of inertia lij and, therefore, real constant d, the 
solution (|6.ip is real. To get real solutions for the other components one just takes real and 
imaginary parts of (|6.2p . (|6.3p . as well as 72*\ 73*\ using the formula e'"* = cos(at) + i sin at. 

In view of (|6.4p . as i ^ ±00 the real vector solution 7^^' starts and ends along the axis 
(1, 0, 0), although for finite t its evolution can be complicated (see Figures 4 and 5 below). 
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6.2 Examples of meromorphic solutions. 

Below we show two examples of independent meromorphic solutions {71,72173} of the 
Poisson equations for p = 1 and p — 3. First, the solutions (|6.ip . (|6.2p were taken, their 
real and imaginary parts were extracted, then the formulae (|6.5p were applied. We write the 
obtained real solutions in terms of hyperbolic and trigonometric functions. 
The simples solutions for p — I are 

7 -(^tanh(t), 

(2) ^ / ( i\ gg^j^(^) cos(;|)tanhft)-(isin(|) rfcos (|) tanhffl + sin(|) ^ 
y \dj ' ^(F + 1 ' y/(F + 1 

7(3) = fsin (-) sech(i) + '^co^(l) c;sin(i) tanh(t) - cos(i) 



dj Vd^TT Vd^TT 

where sech (t) = 2/(e* + e~'). 

As one can check, these vectors form an orthonormal basis, (7^*\ 7''-'-') — Sij, and, as 
t ±00, we have 7^^^ = (±1,0,0) and, respectively, 

(2) /^Q cos{t/d) — dsin{t/d) d cos{t / d) + sin{t / d) 
(•j) cos(i/d) + (isin(i/(i) (icos(t/(i) — sin(t/(i) 

Since d > 0, as t ^ —00 (respectively t — > —00), in the body frame the vectors 7(2)^^^(3) 
perform uniform rotations in the plane (0,1,1) in counterclockwise (resp. clockwise) direction. 

This implies the limit spatial motions of the body are rotations about the same axis, 
with the same angular velocity 1/d and in the same direction, however the rotation angle 
undergoes the phase shift, which is computed to be arccosj^^. An example of the above 
solutions is illustrated in Fig. 4. 
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Figure 4: Independent and orthogonal vector solutions 7'^' and 7'^-* for p = 1 and d = 1/2 in the body frame 



For p = 'i the corresponding solutions are more complicated: 

7^^) = a(^\/d^ + l(^9 + (f - 4:(f sech(i)2^ tanh(t), [Acf sech (t)'^ - 9 (cf + 1)] sech(i), 

d [3{d^ + 1) - sech {tf] sech (t) 
,(2) ^ a .(^_2^/rfr^e*[(d2(3-10e2* + 3e4*)-9(l + e2*)2)cos 

+ 12d(e'** - l)sin 

'3t 



4e' 



3t 



l(^{d'^ - 15) cosh(2t) + 9 + d^) cosh(i) sir 



+ ( (9 - Td^) cosli(2t) + 9 + sinh(i) cos 



3i 



4e 



3t 



3t 

(7d^ - 9) cosh(2t) - 9 - 17d^ ) x cosh(i) sin ( — 



+ d(^{d^ - 15) cosh(2<) - 15 - 7d^^ sinh(i) cos 



27 



and 



v(3) 



+ 12fi(l - e"**) cos 



3t 



4e 



3t 



^((d^ - 15) cosh(2t) + 9 + ^2) cosh(t) ( 



3t 



{7(f - 9) cosh(2i) -9 -(f) sinh(t) sin 



3t 



4e 



3t 



((9 - 7d^) cosh(2i) + 9 + 17d^^ cosh(i) cos 



d( (d^ - 15) cosh(2t) - 15 - Td^ ) sinh(t) sin I ^ 



where 



x/d2TT(d2 + 9) 

These vector solutions also form an orthonormal basis and the corresponding spatial motion 
of the body is similar to that of the previous case p ^ I. An example of these solutions is 
illustrated in Fig. 5. 



g3(t) 





92% o:^^^^ 0.5 g1(,) 
Figure 5: Vector solutions 7^^^ and 7'-^-* for p = 3 and d = 1/2 
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6.3 First integrals. 

As was mentioned in the beginning of Section 2, if = (Pi, ^2,^3) is a single- valued vector 
solution of the Poisson equations and, moreover, the components of P{t) are single-valued 
functions of uji{t),uj2{t), then the system (j2.ip . (|2.2p admits an additional first integral 

F3 = Pi(wi, ^2)71 + ^2(^1,^^2)72 + ^3(^1, '^2)73- (6.6) 

Since this system is homogeneous, then Pi are homogeneous polynomials of a certain degree 
k. 

Now we use the meromorphic solutions 7^^^ obtained in Section 6.1 for odd p to construct 
the corresponding extra integrals. In the simplest case p — 1, comparing the components of 
7^^^ with the solutions (|4.ip for wi, a;2 and recalling the definition of p, d, we easily obtain 

a 

_ Aawijt) _ d 

P2{t) - — ^ - ^T^^2, 

4uj,{t) d^ 

^3(^)=-^ = -^-2, 



which yields the integral 



d d 



'2 



F3 = duJiJi + --p—-jUJ2l2 - ^2 ^ -^ ^273- 

In case of generic odd p the extra integral can be written by means of the special poly- 
nomial solution (|5.14p of the generalized hypergeometric equation (|5.7p with the parameters 
defined in (j5.9p . To do this, it is convenient to define the following homogeneous function 

Q^Qiu;,,u;2) := f[''-^'^'^Ti{z), (6.7) 

where ^ 

F,:=(d^ + l)u:l + ul z:=^, (6.8) 

-Ti 

and ^i{z) is the special polynomial solution (|5.14p of (|5.7p . that is, 

(2j- - 1)!! d^^p^ - 1) ■ ■ ■ (p^ - (2j - 1)2) ^ 

^^'^ f-{ (2j)!! (d2+p2)...((2,^-_l)d2+p2) ^ 

The above definitions imply that if p is odd natural number, then Qi G M[ciJi,ciJ2] is a 
homogeneous polynomial of degree {p — 1). Moreover, (5i(wi, W2) is an even function of wi, 
as well as even function of L02. 

Theorem 6.2. IJ p is odd natural number, then there exists there homogeneous polynomials 
Pi,P2,Pi G M[aji,W2] of the same degree p such that (|6.6p is a polynomial first integral of 
the Poisson equations. Moreover 

Pi - ujiQ, (6.9) 
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and 



P3 
P2 



1 dPi 



p \d'^ + 1 duj 



-LJ2 



dPi 
duo- 



1 dP. 



duJ2_ 

p LUi \d^ + 1 duJi 



UJ2 



duj- 



-UJl 



UJ2 



(6.10) 



Pi 



Proof. As follows from the definition of Pi, the expressions (|6.10p for P3 and P2 are poly- 
nomials. This is clear for P3. Let us show that also P2 G R[Li;i,cj2]. To do this we observe 
that 



dPs 



dQ 



-2pduj2^ pdwiLJ2-„ 2 

ouox OUJl 



+ pd{l + d^)uJi 



- UJl 



dLUlduJ2 



(6.11) 



Since Q is an even function of uji, we have 



dQ 
duo I 



llIiQ, where Q e R[oji,oj2], 



and thus, from (|6.1ip . we have also 

dP. 



— — = W1P3, where P3 G R[a;i, 

OUJl 



Now from (|6.10p it easily follows that P2 G M.[uji,uj2]. 

It remains to show that with Pi, P2 and P3 given by (|6.9p and (|6.10p . the function (|6.6p 
is an integral of the system. Indeed, if it is a first integral, then its time derivative vanishes, 
i.e.. 



= P3 = -E 



37^—7-^^2 - -^t^l ) i^2li + P3^27l - P3W172 + (P2t^l - Pi ^2)73- 

+ 1 OUJl 0UJ2 



The right hand side of the above equation is a linear form in 7^, hence we have the following 
system of three partial differential equations 







aPi 


dPi 


P 




7: UJ2 — 

1 OUJl 


^ — ^1 

0UJ2 


i 




dP2 


dP2 


P 




1 OUJl 


— ^1 

0UJ2 






dP3 


aP3 


P 




— — ^2 — 
1 OUJl 


— ^1 

0UJ2 



P3 = 0, 

2 - P3W1 



(6.12) 



0. 



The form of this system shows that P2 and P3 can be expressed in terms Pi and its partial 
derivatives. The explicit form of these expression is given by (j6.10p . Using them we eliminate 
P2 and P3 in (|6.12p and obtain one partial differential equation for Pi. Next, taking into 
account the fact that Pi is assumed to be a homogeneous polynomial, we obtain an ordinary 
linear equation for Pi. It remains to show that one of its solutions is given by formula (|6.9p . 
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Introduce a new variable to = 102/^1 and define 

Pi{uJi,uj2) = ujlPi{l,uj2/uJi) =: uj''lpi{uj), i = 1,2,3. 

Then, we have 

dPi p-ir n p-i / ■ 1 o o 

— = LJ[ [ppi - Cjp,;J, — UJ{ Pj, « = 1, 2, 3, 

where prime denotes the derivative with respect to uj. This imphes that our system of 
partial differential equations (j6.12p on Pi can be written as the system of ordinary differential 
equations on pi 

~ d{uj'^ +(f + l)p[ + dpujpi + p{(f + l)p3 = 0, 

d(cj2 + d'^ + i)u;p'2 - dpLo^p2 + p{d^ + 1)P3 = 0, (6.13) 
+ d^ + l)ujp'^ - dpoj'^ps, +p{d^ + l)ujpi - p{d^ + l)p2 = 0. 
Resolving this system with respect to P2,P3, we find 

+ 1 + cj2) , d 

^ pid^ + 1) WTTfP'^ 



d^ l + d^ + ujY ^ ^„ 

{1 + d^yp- 

2d^{p - + d^ +L0^) , / d'pil + d' -{p-1)lo'^) 



P2 = ,1 , ^2X2^2 Pi 

Substituting the above expression into the second equation in (j6.13p . we obtain the following 
third order linear equation 

p'l + qip'l + <i2p'i + qspi = 0, (6.15) 

1 3{p - 2)uj 

UJ 1 + rf2 + ' 
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d^LU^-{l + d2)(5p - 4) + (p - l)(3p - 8)^2) ^ (1 ^ ^2)2(1 ^ ^2y 



d2cj2(l+d2+w2)2 
r .2_/^, , ,2n2 , j2n/ i\,,2 , /o , / ^N„N,,4\ 



d2cj(l + rf2+^2)3 
(l+d2)2(d2_(p_l)(l+^2)y 



d'piil + d'Y - 4(1 + d'){p - l)uj' + (3 + (p - A)p)uo'') 



Now, making change of independent variable to 

, ,2 



^2 + ^2 + 1 ~ (e' + e"*)2' 

and using 

dz_2w(z-l)2 d^z _ 2(z - l)2(4z - 1) d^z _ 24tj(z- l)3(2z- 1) 

~ d2 1 ' dZ^- (i2 ^. 1 ' (d2 + 1)2 ' 
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we transform the equation (|6.15p into 



P'l + bip'l + b2p[ + bspi =0, ' = ^ 

2 2 + 3p 
oi = - 



b. 



3 



z 2(z-l)' 

_ (fil + {p^ -6-8p)z + {A + 3p{2+p)-p^)z^)~p'^{z-l) 
^ ~ Ad^z-lfz^ ' 

p'^{-l+p + z-pz) + (f{p^{z - 1) - Ap'^ z + p^ z^ + p{l + (2 -p2(z - l))z)) 



8d2(z- 1)3^2 

Now under the change of the dependent variable 

pi{z) = {z-l)^v(z), (6.16) 

the equation for v(z) becomes exactly the generalized hypergeometric equation (15. 7|) defining 
the function 3F2 with the parameters ()5.9|) . On the other hand, it easy to notice that, up 
to a multiplicative constant, the dehomogenized Pi given by ()6.9p and expressed in terms of 
the variable z coincides with (|6.16p . □ 



6.4 Conclusion 

As follows from the results of Sections 3 and 5, when /13 — and the parameter p us an even 
integer, an interesting situation takes place: all the solutions of the Poisson equations are 
meromorphic, but the equations itself are not solvable in the class of Liouvillian functions. 
In particular, they do not possess extra meromorphic first integrals. In this case it is natural 
to expect that the Poisson equations are reducible to one of the Painleve equations. 

It also should be emphasized that we managed to reduce the Poisson equations to the 
hypergeometric equation under the restriction J13 = 0. The study of solutions of the third 
order equation (j5.2p in the general case, especially of its monodromy group, is an interesting 
open problem. 
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